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Abstract 

We introduce spatiotemporal solitons of the two-dimensional complex Ginzburg- 
Landau equation (2D CCQGLE) with cubic and quintic nonlinearities in 
which asymmetry between space-time variables is included. The 2D CCQ- 
GLE is solved by a powerful Fourier spectral method, i.e., a Fourier spatial 
discretization and an explicit scheme for time differencing. Varying the sys- 
tem's parameters, and using different initial conditions, numerical simula- 
tions reveal 2D solitons in the form of stationary, pulsating and exploding 
solitons which possess very distinctive properties. For certain regions of 
parameters, we have also found stable coherent structures in the form of 
spinning (vortex) solitons which exist as a result of a competition between 
focusing nonlinearities and spreading while propagating through medium. 

Keywords: Ginzburg-Landau, nonlinear waves, vortex solitons, spinning 
solitons, pulsating solitons, erupting solitons, dissipative solitons 



1. Introduction 

The complex cubic-quintic Ginzburg-Landau equation (CCQGLE) is the 
canonical equation governing the weakly nonlinear behavior of dissipative 
systems in a wide variety of disciplines [7]. In fluid mechanics, it is also 
often referred to as the Newell-Whitehead equation after the authors who 
derived it in the context of Benard convection [7J [H] . As such, it is also one 
of the most widely studied nonlinear equations. Many basic properties of the 
equation and its solutions are reviewed in [5J, [6], together with applications 
to a vast variety of phenomena including nonlinear waves, superconductivity, 



Preprint submitted to Elsevier 



February 8, 2013 



Bose-Einstein condensation, liquid crystals and string theory. In particular, 
in nonlinear optics it describes the pulse generation and signal transmission 
through an optical fiber. 

An important element in the long time dynamics of pattern forming sys- 
tems is a class of solutions which we call "coherent structures" or solitons. 
These structures could be a profile of light intensity, temperature, magnetic 
field, etc. A dissipative soliton is localized and exists for an extended pe- 
riod of time. Its parts are experiencing gain/loss of energy or mass with the 
medium. 

In contrast to the regular solitary waves investigated in numerous inte- 
grable and non-integrable systems over the last three decades, these dissipa- 
tive solitons are not stationary in time [2]. Rather, they are spatially confined 
pulse-type structures whose envelopes exhibit complicated temporal dynam- 
ics. Another relevant feature of dissipative systems is that they include en- 
ergy exchange with external sources. Thus, energy and matter can flow into 
the system through its boundaries. As long as the parameters in the system 
stay constant, the structure that could evolve by changing shape exists in- 
definitely in time. The structure disappears when the source is switched off, 
or if the parameters are moved outside of the range of existence of the soli- 
ton solutions. Since these dissipative systems include energy exchange with 
external sources, they are no longer Hamiltonian, and the solitons solutions 
in these systems are also qualitatively different from those in Hamiltonian 
systems. In Hamiltonian systems, soliton solutions appear as a result of bal- 
ance between diffraction (dispersion) and nonlinearity. Diffraction spreads 
a beam while nonlinearity will focus it and make it narrower. The balance 
between the two results in stationary solitary wave solutions, which usually 
form a one parameter family. In dissipative systems with gain and loss, in 
order to have stationary solutions, the gain and loss must be also balanced. 
This additional balance results in solutions which are fixed. Then the shape, 
amplitude and the width are all completely fixed by the parameters of the 
dissipative equation. However, the solitons, when they exist, can again be 
considered as "modes" of dissipative systems just as for non-dissipative ones 

When a system is integrable, its solutions can be constructed using the 
inverse scattering technique [I]. Further generalization to non-integrable sys- 
tems shows that traveling wave solutions can be found and analyzed using 
some basic approaches from nonlinear dynamics. One of the powerful tech- 
niques that can be applied to these systems is that which employs Hamil- 
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tonian equations. For dissipative systems, there are no conserved quantities 
that can be used to solve particular problems. 

Recent perturbative treatments based on expressions about the nonlinear 
Schrodinger equations are generalized to perturbations of the cubic-quintic 
and derivative Schrodinger equations. The cubic Ginzburg-Landau admits 
a selected range of exact soliton solutions [9J. Even the known special solu- 
tions have exotic dynamic responses to external conditions, so that sensitivity 
analysis of the parameters is paramount. These exist when certain relations 
between parameters are satisfied. However, this certainly does not imply that 
the equations are integrable. In reality, general dissipative nonlinear PDEs 
cannot be reduced to linear equations in any known way, therefore an insight 
to the type of solutions that the equation has may be based on numerical 
simulations. To investigate numerically the two dimensional dissipative soli- 
ton solutions of the CCQGLE, and analyze their qualitative behavior, we 
used a powerful Fourier spectral method. 



2. 2D CCQGLE 

The 2D CCQGLE is solved by spectral methods, i.e., a Fourier spatial 
discretization and an explicit scheme for time differencing. In addition, it 
is also crucial to analyze the dependence of both their shape and stability 
on the various parameters of the CCQGLE, viz. the nonlinearity, disper- 
sion, linear and nonlinear gain, loss and spectral filtering parameters of these 
new classes of dissipative solitons, since the solutions experience interest- 
ing bifurcation sequences as the parameters of the CCQGLE are varied. A 
rich variety of evolution behaviors, which include pulsating, exploding and 
spinning dissipative soliton dynamics is revealed. 

In dimensionless form, the 2D CCQGLE with the corresponding cubic- 
quintic terms is a modified NLS and it takes the form jH |2] 

d t A = eA + (h + i Cl )V\A - (6 3 - icz)\A\ 2 A - (b 5 - tc 5 )\A\ A A, (2.1) 

where Vj_ = ^J^( r J~:) is the transverse Laplacian for radially symmetric 
beams. The present study will confine itself to the spatially infinite system 
in two dimensions and will focus primarily on a spatiotemporal behavior of 
dissipative solitons without imposing radial symmetry on the beams, and 
hence the Laplacian V]_ = ^ + depends on both variables x and y. 



The interpretation of the system's parameters in (2.1) depends on the 
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particular field of work. In optics, A(r; t) is the normalized envelope of 
the field, r is the radial coordinate, t is the propagation distance or the 



cavity number. The system's parameters of (2.1) are: e linear loss/gain, 
b±- angular spectral filtering, c\ = 0.5- second-order diffraction coefficient, 
6 3 - nonlinear gain/loss, C3 = 1- nonlinear dispersion, 65- saturation of the 
nonlinear gain/loss, and C5- saturation of the nonlinear refractive index. 

3. Fourier Spectral Method 

The Fourier transform of a function u(x, y) is defined by 

1 poo poo 

T{u)(k x ,k y ) = u{k x ,k y ) = — / e^ x+k ^u{x,y)dxdy, (3.1) 

^ J -co J — oo 

with the corresponding inverse 

POO fOO 

27T 



F-\u){x,y) = u(x,y) = ±- I I e i{k * x+k y y) u(k x , k y ) dk x dk y . (3.2) 



— 00 j —00 



The function u(k x , k y ) can be interpreted as the amplitude density of u 



for wavenumbers k x , k y . 



First, we write the Fourier transform of equation (2.1) as 



A t = [e - (h + ic x ){kl + k 2 y )] A - (63 - tc 3 )\A\ 2 A - (65 - ic 5 )\A\*A, (3.3) 



and rearranging the terms in (3.3), we consider the following ordinary differ- 



ential equation in the Fourier space 

A t = a{k x , k y )A + 0\aJa + jIA^A, (3.4) 
with a(k x , k y ) = e - (bi + ici)(kl + k y l ), /3 = -(b 3 -ic 3 ) and 7 = -(b 5 -ic 5 ). 



In the Fourier space, (3.4) contains the linear part a(k x , k y )A, and the non- 
linear part /3|A| 2 A + 7|A| 4 A 

We solve the initial value problem for the above ODE numerically using 
an explicit scheme i.e., 4 th order Adams-Bashforth for the non-linear part 
and the exact solution for the linear part, which is given by 

A{t) = A{^0)e a{kx ' k y )t . (3.5) 



4 



3.1. Spatial discretization (Discrete Fourier Transform) 
We discretize the spatial domain [— L/2, L/2] x [-L/2, L/2] into n x n uni- 
formly spaced grid points = (xi,yj) with Ax = Ay = L/n and n even. 
Given A(Xij) = A^, i,j = 1,2, ••• , n, we define the 2D Discrete Fourier 
transform (2DFT) of A as 



A kxky = Ax Ay ^ ^ e-^+M^.., /,,.. /,,_ 

i=l j=l 



+ 1, 
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(3.6) 



and its inverse 2DFT as 



A _ 1 
iJ (2tt) 2 



n/2 



n/2 



E E i,j = l,2,...,n. (3.7) 



fc ;c =-n/2+l fc a =-n/2+l 



In (3.6)-(3.7) the wavenumbers fc^ and and the spatial indexes i and j, 
take only integer values. 

3.2. Temporal discretization (4 th order Adams- Bashforth) 

Given t max we discretize the time domain [0, t max ] with equal time steps 
of width At as t n = nAt, n = 0, 1, 2, ■ ■ • , and define A n = A(t n ). Initializing 

A n = we compute the nonlinear terms A3 = T I J 7 (A n ) J 7 x (A n ^ 



and A/5 = J 7 f 7- x {A n ) T'\A-)\ and advance the ODE (|3.4|) in time with 

time step At using the 4-step Adams-Bashforth method. 
The discretized ODE has the following form 



At 
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55/(A n ) - 59/(A n ~ 1 ) + 37 f(A n ' 2 ) - 9f(A n - 3 ) 

with exact solution of the linear part A n+l = A n e a ^ ki " ky ' t , and the nonlinear 
part / defined by f(A) = /3A/i + 7A2. 

3.3. Numerical implementation 

The language used is Fortran90, with gfortran compiler. All the parameters 
needed by the program are written in a datafile so that they can be easily 
modified, without systematic compilation. Among them, most important 
are the parameters of the 2D CCQGLE, which are shown in Table [TJ For 
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2D solitons 


e 61 ci b 3 c 3 b 5 c 5 


stationary 


-0.045 0.04 0.5 -0.21 1 0.03 -0.08 


ring 


-0.1 0.1 0.5 -0.8 1 0.04 -0.02 


pulsating 


-0.045 0.04 0.5 0.37 1 0.05 -0.08 



Table 1: Sets of parameters for 2D solitons 

the ring solitons from §4.4 we have used the parameters of [3], while for the 
stationary and pulsating, the parameters are from [T2] . 

In this paper we are recalculating the results published previously in 
[3j [12], and using this scheme, we hope to find other regions of parameters 
where stable 2D solitons exist, as well as solitons with arbitrary vorticity 
(m). 

An important subroutine is the Fast Fourier Transform (FFT) [TOf ITS] . 
This subroutine uses an algorithm that reduces computational time from 
0(n 2 ) to 0(n log ri) floating point operations. Moreover, it only works for 
periodic functions hence, caution must be made when we select the size of 
the domain for different cases of the soliton solutions. This must be chosen 
in such way that the solitons do not touch the boundaries, otherwise if the 
non-zero part of the solutions propagates through some boundary, it will 
spread throughout the domain from the opposite side. 

The Fortran program ran on the Zeus cluster at Embry-Riddle Aero- 
nautical University. This parallel computer of 256 nodes is rated to over 1 
Tera-flop. Jobs were submitted thanks to Platform Lava installed on the 
cluster. This software manages to distribute all the jobs to different nodes, 
considering priorities. Once the output file are written, we used Tecplot 360 
for visualization. 

The duration of computations depends on the equation parameters and 
mesh size. The most appropriate size to use was 512x512, but other types of 
meshes were considered. Very fast calculations took almost two hours, while 
longer ones took several days. 

4. Numerical Simulations 

4-1. Initial conditions 

Localized structures with vorticity are obtained when initial conditions 
have radial symmetry, and hence we assume that the initial pulse can be 
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reasonably well approximated by a bell-shaped, i.e., 

i. Gaussian shape, 

A(x,y;0) = A e- r2 (4.1) 

ii. Ring shape with rotating phase, 

A(x, y; 0) = A Q r m e~ r ' e imB (4.2) 

where m is the degree of vorticity, A is a real amplitude that should 
generate sufficient power to place the initial condition into a basin of 
attraction of a 2D soliton, and 9 = tan -1 (^f ) is the phase. The widths 
of each of the initial conditions could be either circular (radially sym- 
metric) or elliptic, and are controlled by the inverse widths parameters 
o x and a y , with r = \j (a x x) 2 + (<7 y y) 2 , see Fig. [IJ In this paper we 
have used the vorticity m — 1, but family of solitons of higher vorticity 
have been found [14"] . 




Figure V. Initial shapes of the field. Left: Gaussian shape. Right: Ring 
shape with vorticity m = 1. 

4-2. Parameters 

The existence of stable solutions is strongly linked to the parameters of the 
equation. One of the most important issue is to find the continuous range of 
existence of solitons of the same type in parameters space. Moreover, several 
solutions can coexist for the same set of parameters. These solitons are not 
necessarily stable. Hence, stability is controlled by the parameters of the 
equation and by the choice of the initial conditions. 
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The main parameter of the solution that we monitor in our simulation is 
the energy [12], Q{t) where 



Q(t) 



/oo poo ,b IL 

/ \A(x,y;t)\ 2 dxdy = ^ ^ |A ?J | 2 AxAy (4.3) 
-oo J-oo i=1 =1 



and is a conserved quantity (a dynamical invariant of ( 2.1[ )). For a localized 



solution Q is finite and changes smoothly while the solution stays within the 
region of existence of the soliton. When Q changes abruptly then there is 
a bifurcation and the solution jumps from a branch of solitons that become 
unstable to another branch of stable solitons. Thus, monitoring Q allows 
us to find bifurcations in an easy way. Observing a finite Q also reveals 
the stability of a solution. As soon as the solution becomes unstable, Q 
diverges. In the case that the resulting solution is pulsating, if the solution 
stays localized, the energy will evolve periodically and will remain finite. 
When the optical field spreads out, the energy tends to infinity, and when 
the solution dissipates the energy evolves to zero [12J. 

4-3. Stationary solitons 

The initial shape is a circular 2D Gaussian with Aq = 2.5, and a x = a y = 
1. Steady solitons are the simplest solutions we can find. After a transition 
period of stabilization, the shape of the solitary waves remains the same. 
The energy becomes constant after z ~ 40s. The soliton profile is always a 
Gaussian, and each profile has a plain bell-shape. Finally, energy converges 
to Q(t) ~ 18.6, and stabilizes forever if the system's parameters are not 
changed. 

4-4- Ring vortex (spinning) solitons 

We choose a circular ring type structure with vorticity m = 1, A = 2.5, 
and a x = o~ y = I, and parameters from line 2 of Table [I] By keeping all 
the parameters fixed, we can find different types of ring vortex solitons, by 
varying the initial conditions. This also can be achieved by varying slowly 
one of the equation parameters while the initial condition stays fixed. First, 
we found a radially symmetric ring vortex soliton shown in Fig. |3j Even 
though the structure oscillates before converging to a stable value of the 
energy as it is with the stable soliton, now since we introduced a phase, the 
soliton is spinning. The rotation has a fixed angular velocity and can occur 
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in either direction clockwise or counterclockwise, which is determined by the 
initial condition. 

Next, we introduced for the same set of parameters, a different initial 
condition, viz. an elliptic shape which breaks the symmetry, o~ x = 0.15, 
cr y = 0.85. This new initial condition leads to another kind of soliton with 
an elongated shape and two peaks diametrically opposed at the top of the 
structure. 

The results for both radially symmetric and elliptic solitons are shown in 
Figs. [3] and [4} 

Ordinary bell-shaped solitons with zero-vorticity also exist for the same 
parameters, but with lower energy. If a ring vortex soliton loses its stability, it 
will be transformed into several bell-shaped solitons via multiple bifurcations. 
That is what we have found, by chance, while looking for the pulsating ring 
vortex soliton, see Fig. [5j This class appears to be unstable, may lead to 
chaos, and needs further investigation. 

4-5. Pulsating solitons 

The following results contain two interesting features: pulsations and 
transitions between stable and unstable states. The initial shape is also 
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Figure 3: Circular ring vortex soliton with radially symmetry. Top left: 
Energy. Top right: Ring vortex soliton at t = 20s. Bottom left: Contour 
plot of \A\ 2 . Bottom right: Phase plot that indicates rotation at t = 20s. 



10 




11 




Figure 5: Unstable vortex circular ring. Left: 10 bell-shaped solitons ap- 
peared in results to the destruction of the initial ring soliton. Right: the 
phases of the solitons are not spinning. The parameters of the initial shape 
are: A = 3.0, o x = a y = 0.15. 

Gaussian with A = 5, but slightly asymmetric, a x = 0.8333 and a y = 0.9091. 
Here we have used the parameters from line 3 of Table [T] For the first 200s, 
the evolution of the solution looks similar to the stationary soliton, see Fig|6j 
While the soliton is stabilized at t — 430s, a perturbation or a bifurcation 
appears to make the structure evolve toward a new higher state: the energy 
Q of the soliton is now pulsating. The structure still has radially symmetric 
profile when Q takes its maximum, but elongates alternatively in different 
directions when Q takes its minimum value, see FigjTJ 

Indeed, this state is not stable for this set of parameters and evolves into 
another stable soliton which keeps the elongated shape, as if it was a double 
bell-shaped soliton. The energy still pulsates but with a very little amplitude 
and the soliton rotates around its center. Different types of pulsating solitons, 
with varying energies and shapes, that rotate or not, have been identified, 
but there is no doubt that there are other types exist too, especially when 
two or more solitons coexist in the same plane at the same time. 

4-6. Exploding solitons 

Another class of 2D solitons is the exploding or erupting solitons. The 
evolution starts from a localized initial condition of Gaussian shape with 
parameters A = 3.0, a x = a y = 0.3. After a while its slopes become covered 
with small scale instabilities which seem to move downwards. The soliton 
explodes intermittently, thus resulting in significant bursts of power above the 
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Figure 6: Pulsating soliton with changing state. The energy, on the left 
shows these transitions. Right: bell-shaped soliton t = 200s. 

average, recovering the initial radially cylindrical shape after each explosion 
[I] . Using the same parameters as in the case of the ID case of papers [13j HTJ , 
i.e., e = -0.1,6! = 0.125, c x = 0.5, b 3 = -l,c 3 = 1,6 5 = 0.1, c 5 = -0.6, the 
energy Q for an exploding soliton existing in the 2D case is shown in Fig. 
|8j As it was shown for the ID exploding solitons, this unusual dynamics 
appears as a result of an instability. As we move further from this boundary, 
the explosions become more violent, and, as a result more than one beam can 
be generated in some cases. This gives rise to a very complicated dynamics, 
ending up with the whole numerical grid filled with the solution. 

At first, the shape seems to be stable with constant energy and smooth 
bell-shaped cross section (t ~ 28s). Often, circular waves go from the center 
to the exterior and vanish as if energy were dissipated in the boundaries of 
the solitons (t ~ 32s). Suddenly, and periodically, the soliton "explodes" 
(t « 35s): its shape grows, so do the energy, until it collapses and goes back 
to the original shape (t ~ 36s). This completely chaotic, but well-localized 
structure restores to a perfect solitons shape after the main burst. The pe- 
riodic evolution of the energy, with high bursts of energy can be seen in 
Fig. [8] This process never repeats itself in successive periods, however, it 
always returns to the same shape. Hence, the essential features of the explo- 
sions, both observed experimentally in ID by [TU [3] are: explosions occur 
intermittently, the explosions have similar features but are not identical, and 
explosions happen spontaneously being triggered by perturbations. These 
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Figure 7: During the pulsating state t e [480s, 520s] the soliton elongates 
alternatively in different directions. 
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Figure 8: Exploding soliton: energy is periodic with high bursts almost every 
12s. 

type of solitons need more numerical studies to be understood in their full 
complexity, and in the future we hope to find other parameters of existence 
for the 2D exploding soliton. 

5. Conclusion 

In conclusion, based on the numerical simulations, we have presented 
evidence for the existence of stable stationary 2D solitons in a dissipative 
medium, described by the 2D CCQGLE, in which asymmetry between space- 
time variables is included. For the domains of the parameters that we ex- 
plored, we have also found stable coherent structures in the form of spinning 
solitons. These exist for a certain region of parameters, as a result of a 
competition between focusing nonlinearities and spreading while propagat- 
ing through medium. Beyond the regions of stability, the solitons may lose 
the radial symmetry but remain stable. Varying the parameters, and using 
Gaussian and ring initial conditions, the simulations revealed more distinct 
2D solitons in the form of pulsating, and for the first time two-dimensional 
exploding solitons with very interesting properties. 
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Figure 9: Exploding soliton at t — 28s, the shape is smooth. Then, circu- 
lar waves appear and grow. Finally, soliton explodes around t = 34s and 
collapses after which it restores back to its initial shape. 
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